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SUMMARY 


The results of calculations are presented for the purpose of showing 
the effects of ionization and electronic excitation on the flow properties 
of one-dimensional shock waves in xenon gas. The calculations are based 
on statistical -mechanics theory with the assumption that thermal equilib- 
rium exists at all points in the flow. The calculations are made by both 
including and neglecting the effects of electronic excitation, and the 
results of the calculations are qualitatively compared with experimental 
data. 


INTRODUCTION 


Air ih passing through sufficiently strong shock waves does not 
behave as an ideal fluid because of the effects resulting from molecular 
dissociatioh, electronic excitation, and ionization. Several theoretical 
investigations of such strong shocks have been made by using the methods 
of statistical mechanics (refs. 1 and 2). The experimental investigations 
of strong shocks in air have been hindered by the difficulty in producing 
them. However, in an effort to study these effects, strong shock waves 
have been produced in other gases (refs. 3 to 5 )* 

Resler, Lin, and Kantrowitz (ref. 4) examined the flow of argon in 
shock tubes and obtained both experimental and theoretical results. This 
investigation indicated that no flow effects would be attributed to elec- 
tronic excitation for flow velocities up to a Mach number normal to the 
shock of 20; therefore, the theoretical calculations were simplified by 
neglecting electronic excitation. The work of reference 3 consisted of an 
investigation, by a ballistic method, of strong shock flows about small 
metal cone-cylinder models in xenon, in which higher Mach number flows 
can be obtained than in argon. 

The purpose of this paper is to present the results of theoretical 
calculations of the flow properties of xenon in strong shock waves . 
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Two sets of calculations are made, one including the effects of electronic 
excitation and ionization, designated as excitation calculations, and the 
other including only the effects of ionization, designated as nonexcitation 
calculations. The excitation calculations follow the method described in 
reference 1. The results of the calculations are qualitatively compared 
with experimental results given in reference 3 on the assumption that the 
glow observed in a xenon shock wave is approximately proportional to the 
percent ionization. 


SYMBOLS 


7 

Xe 


adiabatic index, 5/3 

xenon /atomic weight, 131.3 — - — ) 
y mole/ 

Avogadro's number, 6.02JM- x 10 ^ molecules 

mole 


m c 


Boltzmann's constant, 1. 38026 X 10 

-28 

electron mass, 9*10721 X 10 , g 


-16 ergs 


°K -molecule 


energy, ergs/mole 


energy equivalent, 1.60186 x 10 


•12 ergs 

electron volt 


K 

X 

a 

M 


m 

5 


equilibrium constant 

first ionization potential of xenon, 12.078 electron volts 
ionization, percent 
free -stream Mach number 


Mach number normal to shock, 


U 1 

\p Pl V l 


mass 

one-half included angle of cone, deg 
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Z partition function 

h Planck's constant, 6.62377 X 10" 2 7 erg-sec 

p pressure, dynes/cm 2 

Uc/uj ratio of velocity on cone surface to limiting velocity 
0 shock-wave angle, deg 

V volume, — c - m ^ - 

g-mole 

gjj statistical weight 

T temperature, °K 


v n 

J 

R 


term-value difference, cm"l 
total angular momentum of molecule 


universal gas constant, N Q k, 


ergs 

°K-mole 


u velocity, cm/sec 

c velocity of light, 2.99790 x 10 1 ® cm/sec 

Subscripts.: 

I non-ionized or neutral atom 

II ionized atcm 

1 a point in the otherwise undisturbed stream 

2 a point immediately behind shock wave 

e electron 

ex atomic excitation 

i internal part 

ion atomic ionization 
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n the nth term 

t atomic translation, also translational part 

The term mole refers to a mole of gas before entering the shock wave. 

METHODS OF CALCULATION 
Excitation Calculations 

• When there is no electronic excitation or ionization in a monatomic 
gas, the flow properties before and after a normal shock wave can be 
given by the equations of continuity, momentum, and energy and the per- 
fect gas laws. These are as follows: 

continuity 


momentum 


energy 


state 


U 1 V 2 

= u 2 V 1 



(1) 

1! 

u 2 2 

*2 + — 



(2) 

u i 2 = 

e 2 + p 2 v 2 

+ 1 u 2 2 


( 3 ) 

p i v i 

= RT X 



(M 

P 2 V 2 

eh W 

II 



( 5 ) 

and E2 are 

given by 

| RT i 

and ^ RT 2 > 


respectively. If a shock-flow problem is adequately specified, these 
five equations uniquely determine the unknown flow properties. 

For a real monatomic gas, the temperatures encountered in strong 
shock waves induce electronic excitation and ionization which become 
important and make the equation of state (5) and the expression for E 2 

no longer valid. The correction for ionization and electronic excitation 
is made in the following manner: First, the equation of state which is 

true for a constant number Qf moles of particles per unit mass of gas is 
changed to 


P 2 V 2 


■ l 1 + skH 


(6) 
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so as to include the increased number of particles due to ionization. 


where a is the percent ionization in a shock wave when there is no 
ionization of the gas entering the shock. Second, the energy E 2 , which 


reference 1 for a diatomic gas. They are altered here to apply to a 
monatomic gas by dropping terms corresponding to molecular vibration, 
rotation, and dissociation. 

The value of a is determined from the equilibrium constant K 
which is related by Saha's equation to the partition functions and the 
temperature of the gas (ref. l). Thus, 




( 7 ) 


where 


* t 



( 8 ) 



( 9 ) 


E ex = RT 2 2 



( 10 ) 


The energy expressions (eqs. (8) to (10)) are given in section II of 



2 * 


( 11 ) 
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from which there is obtained 


a 



( 12 ) 


In this analysis, the temperatures considered are not high enough 
for double ionization of an appreciable number of atoms. The number of 
doubly ionized atoms is therefore neglected, and a represents the per- 
cent of atoms that are singly ionized. 

The partition functions of the atoms Z which appear in equations (10) 
and (ll) are defined by 


Zj . 


(13) 


where 



( 14 ) 



( 15 ) 


and 


z n - z n t z n i 


(16) 


where 



(17) 



(18) 
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For an electron. 



(19) 


where 



( 20 ) 



( 21 ) 


These equations are standard expressions of statistical mechanics and 
equivalent forms may be found, for example, in reference 1. They are 
based on the assumption of thermal equilibrium. The summations indicated 
by equations (15) and (18) extend over all the electronic energy levels 
of each atom type, that is, the neutral and the ionized atoms. The term- 
value differences measured from the ground states are v T and v" TT 

x n lx n 

and the g's are the statistical weights. For the nth state, g n 
is 2J n + 1 where J is the total angular -momentum quantum number as 
given in reference 6. 

With these revisions, ionization and excitation effects are included 
•in the analysis. The equations to be solved are equations (l) to (4), 

(6), and (ll) and the relations for E 2 given by equations (8) to (21). 

The procedure for the calculations is as follows . 

By using spectroscopically determined term-value differences for 

each atom type, Xe T and Xe TT , up to v T = 97,007.41 cm -1 and 

11 -L n 

Vji n = 157 * 225.10 cm"-*-, respectively (see ref. 6), the internal part 

of the partition functions (eqs. (15) and (l8)) are calculated for values 
of T 2 ranging from 300° K to 3,000° K in intervals of approximately 500° K. 

The variation of these partition functions with temperature is shown in 


tion (10) are determined and are plotted as a function of T 2 (see fig. 2). 


figure 1. From these curves the two RT 2 ^ ~ log terms in equa- 
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The total partition functions from equations (13), (16), and (1'9) are 
then determined for various values of T 2 and are substituted into 

equation (ll) to obtain a as a function of T 2 and V 2 . The ioniza- 
tion a is then plotted as a function of T 2 * For a given value of V 2 , 

the corresponding values of a, RT 0 2 — log Z T , and RT 0 2 log Z TT . 

2 dT I i 2 dT ■ Lli 

are substituted into equations (8) to (10), and the energy E 2 is 

plotted as a function of T 2 * 


In order to satisfy a shock-flow condition, the proper value of E 2 

is determined by use of equations (l) to (3) which relate conditions 
before and after the shock. If u-^, u 2 , and p 2 are eliminated from 

equations (l) to (k) and (6), then 


Eo = 


| FT i"! 


1 + 1 + 


•100 /T x V 2 


— + 3r 


(22) 


Standard conditions are chosen for p-^, V^, and T-^. 'Since a is known 

as a function of T 2 and V 2 , this second expression for E 2 is plotted 

as a function of T 2 for the given value of* V 2 . The two functions 

of E 2 are compared on the same graph for the same value of V 2 . They 

form two intersecting curves which fix the values of E 2 and T 2 for 

one shock condition. From these values the unknown flow conditions u-j_, 

U 2 , and p 2 are determined. Different shock conditions are calculated 

by selecting a different value of V 2 and by following the same procedure. 

The final results are plotted as a function of Mach number of the undis- 
turbed stream. 


Nonexcitation Calculations 


Resler, Lin, and Kantrowitz (ref. k) have shown that electronic 
excitation can be neglected for flows in argon up to M-j^ = 20. Since the 

electronic excitation levels of xenon differ from those of argon, calcu- 
lations are made to show the results when the electronic excitation of 
xenon is neglected. The equations used are equations (l) to (k), (6), and 
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E 2 = E t + E ion 



+ — N 0 Xa 
100 0 


(23) 


to which equation (j) reduces when excitation is neglected. In this 
case the internal part of the partition functions for the atoms in equa- 
tions (15) and (l8) reduces to the first term of the summation series. 
Thus, equations (13) and (l6) become 



(2M 


and 


where 



(25) 


g Io = 1 

and 


as obtained from reference 6. 

Using this result gives Saha f s equation in the following form: 


K = 8 


/2nm e k\ 5//2 V 2 T 2 3 / 2 -Xa/kT 2 


N, 


(26) 


where the mass ratio is approximated as unity. 

In order to determine the flow conditions of a shock wave by this 
method, the percent ionization from equation (12) is plotted as a function 
of Tg for a given value of V 2 and the energy Eg in equation (23) 

is then plotted as a function of Tg. As in the preceding case, values 
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of E 2 and T 2 which satisfy equations (22) and (23) for the given 
value of V 2 are determined graphically. All the properties of the 
flow are then determined as functions of the free-stream Mach number. 


RESULTS OF CALCULATIONS 


The excitation calculations and the nonexcitation calculations 
are made for one -dimensional shock waves in xenon up to M 1 = 15. 

Standard conditions are assumed in the undisturbed flow which are 

p-j^ = 1,013,250 dynes/cm 2 , T = 273.16° K, and = 22,420.8 cm^/mole. 

The values of constants used in the calculations are obtained from ref- 
erences 7 and 8. The results of the calculations, together with the 
appropriate ideal-flow quantities, are presented in graphical form plotted 
against Mach number of the undisturbed flow. Percent ionization, tem- 
perature, specific -volume ratio, pressure ratio, and stream velocity 
behind the shock wave are shown in figures 3, 4, 5, 6, and 7, respectively. 

Generally it is evident that the noticeable effects of ionization 
and excitation begin at approximately M-^ = 8 and increase as the speed 

of the shock increases (fig. 3). When ionization and excitation are 
present, these processes absorb part of the energy of the gas and leave 
less translational energy than for the ideal flow. As an almost direct 
result, the excitation-calculation temperature is less than the ideal- 
flow temperature (fig. 4). The specific -volume ratio and the pressure 
ratio are influenced by two effects, namely, the lowered temperature and 
the change in the number of particles per unit mass of gas (a direct result 
of ionization) . In addition, of course, they are related to each other 
by the equation of state. The net result is that the excitation- 
calculation specific -volume ratio is considerably greater than in the 
ideal-flow case (fig. 5) and that the excitation-calculation pressure 
ratio is somewhat larger than in the ideal-flow case (fig. 6). The 
excitation-calculation stream velocity behind the shock reaches a maxi- 
mum at M 1 * 12, whereas in the ideal-flow case it continues to increase 

(fig. 7). 

A comparison of the results of the excitation calculation and the 
nonexcitation calculation shows that, within the range of the present 
calculations, the effect of excitation is of the order of 10 percent of 
the effect of ionization. When the assumption is made that no energy of 
the gas is absorbed as excitation energy, more energy is available for 
ionization and translation. Thus, the nonexcitation percent ionization 
is larger than the excitation percent ionization (fig. 3) and the non- 
excitation temperature is larger than the excitation temperature (fig. 4). 
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The results of the two methods of calculation for specific -volume ratio, 
pressure ratio, and stream velocity after the shock are in agreement with 
the physically plausible concept that calculations in which excitation 
effects are neglected should give results that are less different from 
the ideal -flow case than calculations in which excitation effects are 
included. The differences between these two calculations increase as 
the flow speed increases . 

Figure 8 shows the turning angle of the stream through two types of 
shock waves which have various inclinations. The excitation-calculation 
results are compared with those of the ideal-flow case. Since the presence 
of excitation and ionization causes the stream velocity to be smaller 
after the shock, the stream turning angle is greater in this case. Larger 
turning angles in this flow indicate that at a given Mach number the same 
angle of shock is obtained by a blunter shaped body. The values of shock 
inclinations about conical-shaped bodies in an axisymmetrical ideal flow 
are given in table I. The shock angles are calculated for a monatomic 
gas where y = 5/3 by the method described in reference 9* 

The effects of electronic excitation can be neglected for flows in 
some gases. Resler, Lin, and Kantrowitz (ref. 4) showed that this situ- 
ation is true for the case of argon since for flow velocities as high 
as M]_ = 20 the amount of electronic excitation energy is approximately 
1.7 percent of the total internal energy of the gas. For a flow in xenon 
this proportion of electronic excitation energy is present in the shock 
wave at M]_ = 11.2 (see fig. 9) and becomes 13.3 percent at = 15. 

Thus, the effect of electronic excitation is much more important in xenon 
than in argon. 


COMPARISON OF THEORETICAL AND EXPERIMENTAL RESULTS 


If the presence of ionization in a shock wave is assumed to be 
characterized by a glowing phenomenon and the intensity of the glow is 
assumed to be approximately proportional to the percent of ionization, 
a qualitative comparison can be made between the results of reference 3 
and theory. For this purpose figures 8 to 10 of reference 3 which show 
flows in xenon about cone-cylinder models are reproduced as figures 10 
to 12. In addition, the flow about a small steel sphere is shown in 
figure 13. Figure 13 was obtained by the apparatus used in reference 3. 
The test results are given in table II. 

The calculations show that ionization is present in a shock wave for 
Mach numbers greater than approximately 8. Figures 10 to 12 indicate 
this occurrence to be approximately true for the actual shock flows. A 
glow is seen in the shock wave of figure 10 where the Mach number normal 
to this shock M^ is 9*0 ana this glow covers the silhouette of the 
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cone part of the cone -cylinder model. For this case, theory indicates 
0.06 percent ionization (see table II). When the Mach number normal to 
the shock is reduced, the intensity of the glow in the shock wave dimin- 
ishes nearly proportionally until it is no longer visible. Figures 11 
and 12 show cases of strong shock flows where the intensity of the glow 
is not observed. In figure 11, the Mach number normal to the shock wave 
is k.9. The theory indicates no appreciable ionization in a shock wave 
of this strength. Likewise, figure 11 shows that there is no glow on the 
silhouette of the cone part of the cone-cylinder model to correspond to 
that in figure 10. The flow about the cone-cylinder model in figure 12 
indicates the approximate Mach number range in which glow begins. The 
model is yawed k° to the free-stream direction. The normal Mach number 
of the shock about this yawed cone varies from 8 to 10. A glow is asso- 
ciated with the greatest shock inclination and no glow is associated 
with the least shock inclination; this result corresponds to theory. As 
would be expected from theory, figure 13 shows a more intense glow than 
figure 10 since 2.9-percent ionization is calculated for the sphere com- 
pared to 0.06 percent for the cone -cylinder model. 


CONCLUDING REMARKS 


The work of Resler, Lin, and Kantrowitz indicates that simplified 
calculations may be used to determine the flow properties of strong shock 
waves in argon because the electronic excitation energy is small compared 
to the total internal energy. These investigators use these simplified 
calculations for flow velocities up to a Mach number normal to shock M^ 

of 20 without large error. In the present investigation it is found by 
•comparison with more exact calculations that the effect of electronic 
excitation is much more important for xenon than argon. The difference 
between excitation and nonexcitation calculations increases with increases 
in flow speed. 

The flow properties which are calculated by considering the effects of 
electronic excitation and ionization are compared with ideal -flow prop- 
erties. Differences occur between these two flows when ionization is 
present. The differences increase as the ionization increases at larger 
flow speeds . The calculations show that the actual temperature and stream 
velocity behind the shock wave are less than in the ideal case, whereas the 
pressure ratio and specific -volume ratio across the shock are greater. 

The real stream velocity behind the shock reaches a maximum value as the 
speed increases and then decreases. 

The results of the calculations are qualitatively compared with 
available experimental results . At shock speeds above approximately 
M-^ = 8 the experimental results show a visible glow in the shock. 
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If the assumption that the intensity of the glow is approximately pro- 
portional to the percent of ionization is used, the experimental results 
appear to he in qualitative agreement with the results of calculations . 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., September 25, 1953. 
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TABLE I 


PROPERTIES OF SYMMETRICAL CONICAL SHOCK FLOW IN AN IDEAL GAS 

[r - 5/5] 


5, 

deg 

«c An 

M 

0 

deg 

Uc/ u Z 

M 

e 

5 

0.650 

.800 

.860 

.870 

.880 

.950 

.970 

.975 

.980 

.990 

1.5283 

2.4048 

3.1900 

3.2202 

3.3898 

5.8263 

8.0133 

9.7652 

10.5961 

18.5376 

40° 58' 
24° 46' 
19° 19' 
18° 25' 
17° 32' 

10° 56' 
8° 52' 
8° 16 ’ 
7 0 41' 
6° 27' 

30 

0.400 

.600 

.650 

.700 

.750 

.785 

.800 

.810 

.814 

.815 

.818 

.820 

1.5700 

2.3327 

2.8481 

3.5369 

4.5525 

6.4434 

8.3343 

11.1318 

13.0189 

13.7233 

17.6227 

21.5499 

66° 8' 
46° 6' 
43° 

41° 1' 
38° 3' 
36° 38’ 
36° 4* 
35° 42' 
35° 34' 
35° 31' 
35 0 26' 
35 0 20' 

10 

.700 
.800 
.850 
.875 
i .922 
.950 
.970 
.972 

.974 

.976 

1.9^98 
2.6190 
3.2402 
3.6968 
5.2171 
7 >371 
12.8992 
14 .4189 
16.4326 
19.8942 

33° 9' 
24° 30' 
20° 43' 
18° 55' 
15° 38' 
13° 41' 
12° 16 ' 
12° 8' 
11° 59' 
11° 50' 

40 

.400 

.500 

.600 

.632 

.650 

.667 

.675 

.676 

2.3018 

2.8061 

4.3203 

5.6989 

7.2907 

11.2260 

17.6869 

20.3139 

65° 4L' 

57° 31' 
51° 41' 
50° 11' 
49° 25' 
48° 43' 
48° 23' 
48° 20' 

20 

.600 

.700 

.800 

.825 

.850 

.860 

.870 

.880 

.890 

.900 

.910 

.916 

.920 

1.7521 
. 2.3300 

3.4521 
3-9399 
4.6710 
5.0613 
5.5861 
6.2865 

7.2219 

8.7875 

12.7897 

16.8710 

28.1718 

42° 58' 
34° 53' 
28° 52' 
27 0 37' 
26° 23' 
25° 56' 
25 0 26' 
25° 3' 
24° 34' 
24° 6' 
23° 52' 
23° 26' 
23° 15' 

45 

.375 
.400 
.500 
.525 
.535 
.550 
.570 
• 575 
.578 
.580 

3.1767 
3.2403 
4.5126 
5-3483 
5.9554 
7.1028 
11.7769 
14 .4478 

18.4953 

22.4356 

68° 56' 
66° 52’ 
59° 55’ 
58° 32’ 
57° 59’ 

570 iL' 

56° 19' 
56° 3' 
55° 55’ 
55° ^9' 
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TABLE II 


EXPERIMENTAL TEST RESULTS OF SMALL CONE -CYLINDER MODELS 
AND SPHERES FIRED INTO XENON GAS 


Figure 

5, 

deg 

Model 

velocity, 

cm/sec 

Tl, 

°K 

Free-stream 
Mach number, 
M 

e, 

deg 

(observed) 

a, percent 
(calculated) 

a 10 

J+5.0 

194,923 

296.4 

n.o 

56 ° 31.0' 

0.06 

a ll 

22.5 

192,116 

297.6 

10.8 

26 ° 45 . 5 ' 

0 

a 12 

45.0 

183,446 

295.8 

10.4 


— 

13 

90.0 

219,456 

294.1 

12.5 


2.9 


a Data obtained from reference 3. 



Figure 1.- Variation of internal partition functions 
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Figure 2.- Variation of electronic excitation energy with temperature . 
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Figure 2.- Concluded. 


NACA TN 5091 


19 



Figure 3>- Percent ionization calculated for shock waves in xenon. 
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Figure 4.- Variation of temperature in shock wave with Mach number normal 

to shock. 
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Specific -volume ratio, 



Figure 5»- Variation of specific-volume ratio with Mach number normal 

to shock. 
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Pressure ratio 



Figure 6.- Variation of pressure ratio with Mach number normal to shock. 
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Normal stream velocity after shock, cm/sec 



Figure 7-- Variation of normal stream velocity after shock with Mach 

number normal to shock. 
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Stream turning angle, deg 
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ro 
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Mach number normal to shock 


Figure 8.- Stream turning angle through a shock wave with and without 

ionization. 
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100, percent 



Figure 9 *- Percent of excitation energy in a shock wave 



L- 81272 

Figure 10.- A shadowgraph of the flow phenomenon about a 90° cone- cylinder 

model traveling at M = 11.0. 
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l-81273 

Figure 11.— A schlieren photograph of the flow phenomenon about a 45^ cone- 

cylinder model traveling at M = 10.8. 
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L-8127U 

Figure 12.- A shadowgraph of the phenomenon about a yawed 90° cone-cylinder 
model traveling at M = 10.4. Angle of yaw, 4°. 
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IBMf 


M 



l - 802 6k 

Figure 13.- Photograph of the glowing area about a steel sphere traveling 
at M = 12.5. Approximate contour of model is added to the photograph. 
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